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We study the effect of a missing spin in a one dimensional S = 1/2 antiferromagnet witfi nearest 
neigfibour Heisenberg exchange J and six-spin coupling Q = 4(jJ using Quantum Monte-Carlo 
(QMC) and bosonization techniques. For q < qc f^ 0.04, the system is in a quasi-long range 
ordered power-law antiferromagnetic phase, which gives way to a valence-bond solid state that 
spontaneously breaks lattice translation symmetry for q > qc. We study the ground state spin 
texture $(r) = {Gf\S'{r)\Gf} in the the Stot = 1/2 ground state \Gfrangle of the system with 
a missing spin, focusing on the alternating part Nz{r). We find that our QMC results for N^ at 
q — qc take on the scaling form expected from bosonization considerations, but violate scaling for 
q < qc- Within the bosonization approach, such violations of scaling arise from the presence of 
a marginally irrelevant sine-Gordon interaction, whose effects we calculate using renormalization 
group (RG) improved perturbation theory. Our field-theoretical predictions are found to agree well 
with the QMC data for q < qc. 

PACS numbers: 75.10.Jm 05.30.Jp 71.27.+a 






I 

O 

o 



> 

en 



o 



X 



I. INTRODUCTION 

The one-dimensional 5=1/2 Heisenberg antiferro- 
magnetic spin chain, with nearest-neighbour exchange 
couplings J is perhaps the simplest important model spin 
system in quantum magnetism. It has not only proved 
useful as a theoretical model for the magnetic properties 
of several Mott insulating materialair— , but has also been 
the subject of many theoretical advances such as Bethe's 
original 'Bethe Ansatz' solution of this quantum many- 
body problem and later field-theoretical treatments that 
applied bosonization techniques to map the system to a 
1-1-1 dimensional bosonic field theory with a so-called 
'sine-Gordon' action, made up of a scale invariant free- 
field part perturbed by a non-linear cosine interactioni^ 
In addition, the renormalization group (RG) analysis of 
the cosine interaction that perturbs the scale-invariant 
free-field action is a paradigmatic example of the treat- 
ment of 'marginally irrelevant' interactions in the neigh- 
bourhood of a well-characterized and tractable scale in- 
variant RG fixed point.—"— 

Such marginally irrelevant interactions can give rise to 
violations of scaling predictions at critical points due to 
the presence of logarithmic corrections that multiply the 
scaling answer. A well known example is the 0{N) criti- 
cal point in four space-time dimensions.— In some other 
cases, such marginally irrelevant interactions give rise to 
additive corrections to scaling, which vanish logarithmi- 
cally slowly. The one dimensional Heisenberg chain dis- 
plays both kinds of effects. For instance, gaps in the finite 
size spectra of the spin-half chain are known to have ad- 
ditive logarithmic corrections that do not affect the lead- 
ing behaviour—, while the temperature dependence of the 
NMR relaxation rate 1/Ti violates scaling expectations 
due to the presence of an additional logarithmic factor in 
its temperature dependence- 
Similar logarithmic violations of scaling, arising from 



multiplicative logarithmic factors that multiply scaling 
predictions, have been argued to exis t^'^d^ in a much 
less well-understood case of a two dimensional S = 1/2 
square lattice Heisenberg antiferromagnet on the verge 
of a continuous quantum phase transitio n-'^^d^ between 
the usual Neel ordered antiferromagnetic ground state 
and a spontaneously dimerized non-magnetic state with 
valcncc-bond order. The underlying critical non-compact 
CPi (NCCPi) field theory that has been proposed^ as 
the continuum description of this transition is not as well 
understood from a RG standpoint, and since the numer- 
ics themselves are also more challenging, there have been 
some differences in the interpretation of these results.— ^^^ 

In our own recent workyi^ we have used extensive nu- 
merical computations to establish the presence of appar- 
ently logarithmic scaling violations in the impurity spin 
texture induced by a missing-spin defect at such a quan- 
tum critical point when the system has the usual 5*1/(2) 
symmetry of spin rotations, and ascribed this effect to 
the presence of a yet-to-be-identified marginal operator 
at the putative NCCP^ critical fixed point. In contrast, 
the corresponding spin texture in a system at an analo- 
gous critical point with enlarged SU{3) symmetrji^i was 
found to obey scaling predictions without any logarith- 
mic violations?^, suggesting that the underlying NCCP^ 
critical point describing this SU{3) transition is free 
of such marginal operators. However, parallel work of 
Kauli^ argues that such marginal operators would typi- 
cally not lead to violations of scaling, and finds an alter- 
native scenario more likely. In this alternative scenario, 
both the SU{2) and SU{3) transitions are described by 
fixed points with a leading irrelevant operator with small 
scaling dimension, and the violations of scaling arise from 
the fact that the quantity being studied depends non- 
analytically on this leading irrelevant operator. 

Here, we try and understand the origins of such mul- 
tiplicative logarithmic corrections to impurity spin tex- 
tures by using the one dimensional Heisenberg antifer- 



romagnct as an example. On the analytical side, we 
work within the bosonization framework and use renor- 
malization group (RG) improved perturbation theory to 
obtain predictions for the alternating part of the spin 
texture in this example. These predictions are compared 
with Quantum Monte-Carlo (QMC) results for a one- 
dimensional chain with nearest neighbour Hcisenberg ex- 
change J and six-spin coupling Q — AqJ. The Hamilto- 
nian for this ' JQ3 model' is: 



H 



N 
i=0 



i,i+l 



QY.p^' 



-lP^ 



i,i+l-f^J+2,i+3-f^i+4,i+5 



P,- 



(1) 



where Py = {j — Si ■ Sj) is the projector to the singlet 
state of the two spin-half variables at sites i and j, both 
J and Q are assumed positive, and we impose periodic 
boundary conditions by placing the system on a ring so 
that site N + 1 + k is identified with site k (the total 
number of spins iV -|- 1 is taken even). 

From our QMC results, obtained using the singlet sec- 
tor valence-bond projection method^^, we find that the Q 
term drives a transition to a valence-bond solid phase at 
qc ~ 0.04, so that the system is power-law Neel ordered 
for q < Qc, and VBS ordered for q > qc- Unlike the more 
well-studied case in which such a transition is driven by 
next-nearest neighbour Heisenberg antiferromagnetic ex- 
change couplings, the present JQ3 model does not have a 
sign problem in standard non-zero temperature QMC cal- 
culations (as well as in the ground state projector QMC 
approach), and can therefore be studied at larger length 
scales and greater precision. 

In order to explore the effects of vacancy defects, we 
remove the spin at site and delete all interactions that 
involve this spin from our Hamiltonian. Since N is odd, 
the ground state of the chain with a missing spin is a dou- 
blet with Stot = 1/2. We focus on |G^), the S^^^ = 1/2 
component of this doublet, and compute the spin tex- 
ture <J>(r) = (S'^(r))^ in this ground state for various 
values of q. This is done using a recently developed 
modification^^ of the singlet-sector projector Quantum 
Monte Carlo (QMC) techniquci^. This spin texture can 
be decomposed as $(r) = $„(r) -I- {—l)''/°'Nz{r), where 
alternating part Nz{r) and a uniform part $u(r) are 
obtained from our numerical data by a suitable coarse- 
graining procedure. 

These numerical results for Nz (r) are compared to field 
theoretical calculations within the bosonization frame- 
work, keeping careful track of the effects of the marginal 
cosine interaction term using one loop RG improved 
perturbation theory. Our basic conclusion is that this 
marginal cosine interaction does indeed lead to logarith- 
mic violations of scaling by introducing logarithmic cor- 
rections that multiply the scaling predictions for N^ in 
the power-law Neel phase. Comparing these analytical 
predictions with our numerical results for q < qc, we find 
good agreement with the data, with the strength of the 
log corrections being larger for q further away from the 
critical point, and vanishing for q = qc-, as predicted by 
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FIG. 1: Spin-spin correlation function at distance N/2 
C{N/2) = {Si ■ Si+n) in the ground state of a periodic chain 
with A'^ spins, multiplied by A*' and plotted against q to serve 
as diagnostic of the quantum phase transition from power- 
law Neel order to valence-bond solid order, as discussed in 
the text. 



the bosonization approach. 

The rest of this article is organized as follows: In Sec- 
tion ini we first summarize our approach to the analyti- 
cal calculation of the ground state spin texture induced 
by a missing spin, give our final predictions for the na- 
ture of the logarithmic violations of scaling, and discuss 
them from a somewhat more general RG standpoint. In 
Section Hill we describe our projector QMC studies and 
compare the numerical data for N^ with our analytical 
predictions to establish our main results. We conclude 
with a very brief discussion regarding the connection be- 
tween our results and earlier work on the effect of va- 
cancies on the NMR Knight shift and the spin structure 
factor. 



II. BOSONIZATION CALCULATION OF 
GROUND STATE SPIN TEXTURE 

A. Preliminaries 

As is well-known, we may model our one dimensional 
magnet by the continuum effective Hamiltonian^ 



where the free field part Hq is written as 



Ho = -j^ d. 



dx 



dx 



and the interaction term Hi reads 

"^ /2(/.(a;) 



Hi = 5- / dxcos 



'0 "'0 



V P 



(2) 



(3) 



(4) 



1.4 
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FIG. 2: L^/^N^r) plotted versus r/L (where L = TV + 1 for 
chains with A*' = 959 and A'^ = 479 spins and open bound- 
ary conditions) and compared with the scaling prediction Fq 
for q — 0.04, the approximate location of the quantum criti- 
cal point separating the power-law Neel phase from the VBS 
ordered phase in the one dimensional JQ3 model. Note the 
data at both sizes fits essentially perfectly to the scaling pre- 
diction with the same prefactor Csc- Also note that the best 
two-parameter fit corresponding to our RG improved pertur- 
bation theory result also gives |eo| = 0, and thus coincides 
with the scaling answer. 



here rp is an ultraviolet regulator defined precisely later 
and 



27ri?2 



1 - Treo . 



(5) 



The last constraint that relates R to the bare cou- 
pling constant eo at scale tq arises from the SU{2) spin 
invariance of the underlying microscopic theoryj^ The 
well-known Kosterlitz-Thouless renormalization group 
theory^! applied to the present SU(2) symmetric case 
yields the flow equation 



de 



dlogL) 



PML)) 



(6) 



with the one loop expression for the beta function being 
given by^ 

/3,(e(i:)) = 2^6^(2.) -i(27r)VW- (7) 

This equation can be solved to obtain the running cou- 
pling constant e{L) at scale L as^ 



(8) 
Note that eo is negative in the power-law ordered antifer- 
romagnetic phase in the present sign convention. 

Within this bosonized formulation, the operator S^ (r) 
at site r — ja is represented as2^ 



S'ir) 



a d(j) .4 , _,,r . 
zttR dr y/ro 



R 



(9) 



Here, the coefficient of the uniform part is fixed by SU{2) 
invariance while the coefficient of the alternating part is 
sensitive to microscopic details: A = \/ac where a is 
lattice spacing of lattice model and c is a pure number 
that depends on the microscopic Hamiltonian. 

Finally, we also recall that the 1-point function S = 

{-^ sin ( ^j^ ) )t of tlic operator -^ sin ( ^j^ j can be 
thought of as a function of L and the running coupling 
e(L) for fixed bare coupling eo and fixed r/L. Thought of 
in this way, it obeys the Callan-Symanzik type equation^ 



/ d 



\ dlnL 



oe 



-lie) 



5(L,e(L)|eo,-)=0. 



(10) 



with the anomalous dimension having the expansion 



7(^) = U (J) <L) 



(11) 



in terms of the running coupling e. As is well-known, 
this can be solved to leading order in e(L) to give the 
following scaling law for S 



S = 



Fq ( eo 



VI \<V) 



(l + e(i)i?), 



(12) 



where -Fq 



and R I 



arc some functions of the ra- 



tio -f and the key point about this formal expression for 
S is that all dependence on the ultraviolet regulator r^ 
has been traded in for a dependence on e(L), the run- 
ning coupling at scale L for a fiow that starts with bare 
coupling eo at scale r^. 



B. Overview 

With these preliminaries out of the way, we now outline 
the strategy used below to calculate the alternating part 
of (^^(r))^. The basic idea is to begin by calculating 
the result for this alternating part using the bosonized 



part of the alternating spin density and bare perturbation 
theory to first order in cq for a finite system of length 
L. As we shall see below, this bare perturbation theory 
result will turn out to depend logarithmically on the value 
of the ultraviolet eutoff vq via a logarithmic ultraviolet 
divergence arising from a first order perturbation theory 
contribution proportional to eq log -r- ■ This logarithmic 
divergence makes bare perturbation theory suspect, since 
a notionally small ©(eo) correction turns out to have a 
logarithmically diverging coefficient. 

To extract useful information from the bare pertur- 
bation theory, it is therefore necessary to appeal to the 
Callan-Symanzik equation for the one point function S, 
and use the fact that S is expected to have the general 
form 



Fn 



eo 



^^tII^Z)) (i + ^(^)^)' 



(13) 



as noted earlier. In order to make contact with our bare 
perturbation theory result, we expand this renormaliza- 
tion group prediction to first order in the bare coupling 
constant: 



S 



mi 



TT L 

0^0 log — 
2 ro 



1- Jeolog- + eoi?(7) 
2 To L 



l + eoR{ 



L' 



(14) 



By comparing with the result of our first order perturba- 
tion theory in cq, it becomes possible to fix the functions 
Fq and R. This strategy gives us the one-loop RG im- 
proved result for the alternating part of (>5'^(r))-^ 



N,{r)=c^/^ 



Fo 



eo 



VL \e{L) 



(l + e(L)i?), 



with 



and 



r 

Foi-^) 



■K sin Or 



(15) 



(16) 



T TT 



with er = ^. 




)COt 



(17) 



In order to cast this expression into an explicitly useful 
form for comparison with numerical results on a chain of 
N sites with lattice spacing a, we rewrite the prefactor 



e(L) 



l + 2^|eo|<ilog(- 



^^°^0°^(^ 



1/4 



(18) 



express e(L) as 
e(L) = - 



eo 



1 



2^hl{log(;^) + ilog(log(A))}' 

(19) 
choose the short-distance cutoff as tq = a, and set the 
length L to L = {N + l)a (see subsection III CI below). 
Eqns P^ . (|16p . ([T7|) with these inputs constitutes a the- 
oretical prediction with two free parameters (the overall 
amplitude c, and the bare coupling eo at the lattice scale), 
and we find below that this provides an extremely good 
two-parameter fit of our numerical data in the power-law 
ordered antiferromagnetic phase of the one dimensional 
JQs model. In addition, the spin texture at q = qc, the 
critical end-point of this power-law ordered Neel phase, 
fits extremely well to the scaling function _Fo, to which 
the more general prediction reduces when eq = 0. 

What do these results tell us about the possible origins 
of such multiplicative logarithmic corrections to spin tex- 
tures at other critical points? To explore this, let us con- 
sider the same calculation of the spin texture, but at a 
different critical point with an irrelevant coupling g with 
small scaling dimension a. In other words, we assume 
that P{g) = —ctg + . . . with a small and positive, and 
7(5) = (50 + ^19 + . . . . In this case, the Callan-Symanzik 
equation would predict that N^ satisfy the scaling law 



N^ (r) = exp 



giL) 



7(g) 

Pig) 



dg]F[-,g{L)] (20) 



for some function F (that needs a more detailed analysis 
to determine) . Using the postulated form of the /? and 7 
functions, one can therefore conclude 



N.ir-) = ^FfL^go/L- 



(21) 



Thus, if the critical point in question has no marginal 
operators, the spin texture will quite generally obey scal- 
ing as long as the scaling function F(x, y) does not di- 
verge as 2/ ^ 0. Conversely, if the critical point in ques- 
tion has a marginal operator, scaling will always be vi- 
olated by multiplicative logarithmic factors even if the 
scaling function F(x, y) is perfectly analytic and well- 
defined in the y — > limit. Indeed, in this marginal case, 
the only way of evading a multiplicative logarithmic cor- 
rection would be to "arrange" for the y — ?► limit of the 
scaling function F{x, y) to have exactly the "right" kind 
of singularity needed to cancel the effects of the multi- 
plicative logarithmic correction coming from the expo- 
nential prefactor. One may therefore conclude that un- 
less the scaling function has a particularly "fine-tuned" 
form, scaling predictions for Nz will be generically vio- 
lated by multiplicative logarithmic corrections in the pres- 
ence of a marginal operator. Conversely, irrelevant oper- 
ators can lead to violations of scaling only if the scaling 
function has a divergence as this operator renormalizes 
to zero. 



C. Details 

When a missing-spin defect is introduced into a peri- 
odic spin chain of A^-|- 1 sites, it converts the system into 
a spin chain of N spins obeying open boundary condi- 
tions. These open boundary conditions can be modeled 
by refering back to the original periodic system and re- 
quiring that the spin density is constrained to go to zero 
at the missing site. As is well knownj^i^ this boundary 
condition can be incorporated by expanding the bosonic 
field (t> in terms of bosonic normal modes as follows: 



Hr)^nR^fr^±':^^^mtl±^ (22) 



(r) = cj)o + iJ2 



n=l 

cos 



K^)(«n-4) 



ri=l 

Here, the non zero bosonic commutation relations are 

0(0), go 

rom 
form 



i,\arn,a 



6mm Hq can be written (apart 



from an (infinite) constant |- J2^^i 'T') ^^ the canonical 



Hn 



2 L 



EfumT\ ^. 



(23) 



Thus the ground state jGo) of the unperturbed Hamil- 
tonian is the vacuum for all the a„, and an eigenstate 
of the zero mode go- Indeed, go|Go) = ■kR\Gq) for the 
Stot = 1/2, S'fo( = 1/2 ground state that we wish to 
model (more generally \Gq) is an eigenstate of q^ with 
eigenvalue 2Tr RS^^^). 



Now, the ground state corrected to first order in eq can 
be written formally as 



IG)^IGo)-E(i^)l^)- (24) 



k^Go 



Go 



Here k = {Nn} with n = l,2...oo and Nn being the 
number of bosons in mode n. For an arbitrary excited 
state, we have the unperturbed energy 



i^°({A^4) = ?f +E-"f^" + J 



(25) 



with LOn = iiji, which gives us the following expression 



for the energy denominators: 

E"i{N„})-E"g=J2^nNn. 



(26) 



As a result, our formal expression for the ground state 
corrected to first order in eo now reads 



|G) = |{A^„ = 0}) + 



ueo 



E 



({JV„}|/;cos(if^)|{0}) 



° {W..}#{0} 



^E„T^n 




This gives the following formal expression for the one 
point function: 



^-({0}|^sm(^ 
Vro \ R 



J 



im> + 3 E 



({0}|-^sin(^) |{iV„})({iV„}|/>xcos(H^) |{0}) 



eo 



E 



{JV„}#{0} 



{JVn}#{0} 



({JV4|-l=sin(^)|{0})({0}|/„"dxcos 



2't>{x) 



V 2HI AT 



|{^4> 



(28) 



where wc can set R = l/v27r in the contributions that 
arise from the C(eo) corrections to |Go), as long as we are 
careful to use the full expression R = (27r — 27r^eo)~^/^ « 
{l+n€o/2)/\/2TT when evaluating the first "unperturbed" 
term in order to obtain the latter correct to ©(eo)- To 
evaluate the matrix elements and expectation values, it 
is useful to write the state IJA'n}) in "coordinate" repre- 



sentation as 



{{yn}\{Nn}) 



where the coordinates y„ 




V2 



'HN^iVn, 



(29) 
arc conjugate to "mo- 

th 



menta" 7r„ = °" ^" and H,n{x) is the to*'' Hermite poly- 
nomial of X. The expectation values in our formal pertur- 
bative expression above can now be evaluated in closed 
form using this coordinate representation to obtain the 



following compact integral representation of S 



S {L,-,en 



-eo 



-eo 



TTSmf 



1 



2L 4sin( 
1 



TTSm 



2L 4sin( 



nro 



TT sin 9r f TTEq 

2L \^" ~2~ °^2Lsm9r 



°° r , , ^Sm(0r- 2(1)) /C0s(9r -(!))- COs(er + 

/ dsd(j) 5 — 

Jo sin V coshs - cos(t'r - (p) 



"'0 



dsdcf) 



sm 



ds 



cos 



cos{9r — (j)) 



cosh s — cos( 



(30) 



r 



Here, Or = nr/L, and we have regulated mode sums 
Sm=i 9m over the harmonic oscillator modes by replac- 
ing them with ^^^j^ f/mexp(— TTTTiro/i) whenever nec- 



essary. It is now possible to do the s integrals in closed 
form to obtain the following integral representation for 
S: 



S 



i^^h-) 



2 \2LJ 
2 \2L 



TTSmfc 
2L 



2 °^2isin6l. 



2 sin (h sin dr 
d0 ^ . ; . ^ sin (20 -H 61,-) 



TT — Ur 

1 <.ir 



VSuiS^sin 
2 sin sin ^ 



Vsin Or sin 



" 2 \2l) 



■sin(20-6'r)- 



sin (tt — ((/) — 0r)) 



r 



sin {{(j) + 0,,) - 7r) Jq v^sm^sin 



d0 



2 sin (h sin 0r 



v^suT^sin 
2 sin d) sin 0^ 



sin (2(/) + 61,.) - 



■ sin {2(j) - Or) 



sin (tt — (0 -t- 0r)) 

(0 - 6>^) -K TT 

sin((0 — Or) + tt) 



(31) 



This integral representation is again regulated with the 
short distance cut-off tq by requiring that the integrals 
are to be done by excluding the region [Or — tttq/L, Or + 
■nro/L] from the integration range. Somewhat remark- 



ably, it is possible to obtain explicit expressions for all 
integrals sensitive to this ultraviolet cutoff, and thereby 
reduce this integral representation to the following com- 
pact and simple form: 



J 



S \L,—,€q 



2L 



Treo L 

-TT log — 

2 To 



Treo 




(32) 



Comparing with the general expectation from our RG 
analysis (equation [T4| , we therefore obtain 



r TTsmt 



and 



„, r ^ IT , 27r 
L 2 sm Or 

as already advertised in Section IIIBI 



(33) 




III. NUMERICAL COMPUTATIONS 



Our numerical work on chains with an odd number 
of sites relies crucially on the spin-half sector generaliza- 
tion^^ of the valence-bond projector QMC algorithmi^ 
In our approach, the Stot = 1/2 sector of the Hilbert 
space of an odd number of S* = 1/2 moments, to which 
the ground state belongs, is spanned by a bipartite 
valence-bond cover which leaves one spin 'free'. Roughly 
speaking, the ground state spin texture $(r) = {S^{r))^ 
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FIG. 3: L^''^Nz{r) plotted versus r/L in the power-law or- 
dered Neel phase at g = 0.02 (where L = N + 1 for chains 
with N — 959 and A'^ — 479 spins and open boundary condi- 
tions) and compared with the scaling prediction with a com- 
mon best fit prefactor Csc- Note that the deviation of the 
data from the scaling prediction cannot be simply ascribed to 
an overall multiplicative factor that grows with A'^, since the 
shape of the curves is slightly different. Data at both sizes 
is also fit to the best two-parameter fit corresponding to our 
RG improved perturbation theory result, and the agreement 
is seen to be excellent for the best fit values of crg and |eo 
listed in the legend. 



FIG. 4: L^''^Nz{r) plotted versus r/L in the power-law or- 
dered Neel phase at g = 0.0 (where L = N + 1 for chains 
with A'' = 959 and A*' = 479 spins and open boundary condi- 
tions) and compared with the scaling prediction with a com- 
mon best-fit prefactor Cac- Note that the deviation of the 
data from the scaling prediction cannot be simply ascribed to 
an overall multiplicative factor that grows with N , since the 
shape of the curves is slightly different. Data at both sizes 
is also fit to the best two-parameter fit corresponding to our 
RG improved perturbation theory result, and the agreement 
is seen to be quite reasonable, but not perfect, for the best fit 
values of ceq s.nA |e()| listed in the legend. 



is then obtained directly in our method by keeping track 
of the probability for the free spin to be at various sites 
r (see Ref. [2j for details). 

This method has also been used in computations 
of ground state spin textures at 'deconfined' crit- 
ical points in two dimensional SU{2) and SU{i) 
antiferromagnetsi^i^S, as well as in very recent paral- 
lel work on developing a diagnostic for the presence of 
sharply-defined spinon excitations^l in antiferromagnets. 

We consider pure systems with periodic boundary con- 
ditions and total number of sites ranging from A^ = 60 
to A'^ = 960, as well as the corresponding open spin 
chains obtained by removing one site from the pure sys- 
tem. Our projection power is chosen to scale as 4A^'^ 



to ensure convergence to the ground state. We perform 
> 10^ equilibration steps followed by > 10^ Monte Carlo 
measurements to ensure that statistical errors are under 
control. In systems with a vacancy, we measure $(r) in 
the manner outlined above, and coarse-grain over pairs of 
successive sites to obtain our numerical results for the al- 
ternating part Nz{r), which is to be thought of as living 
on bond-centers in this coarse-graining procedure. We 
have checked that our conclusions are not sensitive to 
the precise coarse-graining procedure used, although non- 
universal details, such as the overall amplitude of Nz{r), 
do change. 

For the corresponding pure systems, we employ the 
singlet sector valence bond projection QMC technique^^. 



and calculate the ground state spin-spin correlation func- 
tion C{j) = (-S'(O) • S{j)) for two sites separated by j — 1 
intervening sites (j < N/2, where N is the total num- 
ber of spins). To begin with, we scan the six-spin cou- 
pling q = Q/AJ and study the q and TV dependence of 
NC{N/2) as a convenient diagnostic that distinguishes 
the power-law Neel ordered phase at small q from the 
spontaneously dimerized valence bond solid (VBS) or- 
dered phase that is stabilized for large q. In the power- 
law Neel phase, NC{N/2) grows (logarithmically) slowly 
with N, while in the VBS phase, it fall off rapidly with 
increasing N. Precisely at the critical point separating 
these two phases, we thus expect a crossing point for 
NC{N/2) plotted against q for various values of N. This 
is precisely what is seen in our data shown in Fig[TJ From 
our data, we estimate that the critical point separating 
these two phases is located at qc ~ 0.04 with an error of 
approximately 0.005 estimated by extrapolating for the 
position of the crossing point (this estimate is consistent 
with the critical point found in Ref . \2T\ \ . 

With this in hand, we compute the ground state spin 
texture in the corresponding chains with one site removed 
for several q < qc for a range of system sizes. The 
alternating part of the computed spin texture is then 
compared with the scaling predictions obtained by set- 
ting eg = 0, as well as with our RG improved pertur- 
bation theory predictions. The former represents a one- 
parameter fit of the data, with the overall amplitude c 
being the only free parameter, while the latter should be 
thought of as a two parameter fit, with the bare value 
eg of the sine-Gordon coupling being the second fitting 
parameter. 

In Fig ^ , we first display our data for the alternating 
part of the spin texture and compare it with the scaling 
prediction at the putative critical point q = qc for two 
of our largest system sizes. As is clear from these two 
figures, the scaling prediction fits extremely well to all 
the data at both sizes. Furthermore, a two-parameter fit 
using the RG-improved perturbation theory result yields 
a best-fit value of eo indistinguishable from eq = 0. This 
confirms our identification of the critical point, since we 
expect that the bare coefficient of the marginally irrel- 
evant cosine interaction is zero at this quantum phase 
transition. 

This excellent fit to the scaling prediction should be 
contrasted with the results shown in Figs (l3]),(|4|), which 
show numerical results at two representative points in the 
power-law Neel phase compared with the one-parameter 
fit obtained from the scaling prediction. As is clear from 
these figures, the scaling prediction simply cannot pro- 
vide a satisfactory account of the data for q < q^ with the 



discrepancy being more pronounced for smaller g, that 
is, further away from the critical point. Furthermore, 
the observed deviations from scaling cannot be simply 
ascribed to an overall N dependent prefactor that grows 
with system size, since the shapes of the curves are them- 
selves slightly different from the scaling prediction. 

In the same figures, we also show the best two- 
parameter fit obtained by using our RG improved per- 
turbation theory result. Two points are worth noting re- 
garding these two parameter fits: Firstly, the best-fit val- 
ues of I eo I increase as one goes further away from q = qc, 
consistent with the expectation that the bare coefficient 
of the cosine interaction vanishes as q approaches qc- Sec- 
ond, the RG improved perturbation theory provides a 
much better fit at g = 0.02 than at the Heisenberg point 
q = — again this is consistent with our expectations, 
since our calculation is perturbative in the renormalized 
coupling e(L), and is therefore expected to provide a bet- 
ter approximation when the bare value of |eo| is smaller 
to begin with. 



IV. DISCUSSION 

We conclude by clarifying the relationship of our calcu- 
lations with earlier calculations of the effect of vacancies 
on2^"— on spin chains. These have typically focused on 
the low- field NMR Knight shift and relaxation rate 1/Ti 
in the presence of vacancies, or the impurity contribution 
to the zero-field spin structure factor and linear suscepti- 
bility of such chains. All these experimental observables 
are obtained from the knowledge of the zero field static 
and equal time spin correlations of the system at finite 
temperature, which has been the main focus of this body 
of work. In contrast, our results focus on local spin tex- 
ture induced by the presence of vacancies at T ~ 0, 
which is a quite different observable connected with the 
impurity contribution to the local susceptibility in the 
high-field regime in which the external field dominates 
over the thermal fluctuations. 
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